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A method to estimate the rehabihty of a perturbative expansion of the stochastic inflationary 
Langevin equation is presented and discussed. The method is apphed to various inflationary sce- 
narios, as large field, small field and running mass models. It is demonstrated that the perturbative 
approach is more reliable than could be naively suspected and, in general, only breaks down at the 
very end of inflation. 
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I. INTRODUCTION 

The stochastic approach to inflation Q, |^ H, Q |^ is an 

efficient method to study how quantum effects can influ- 
ence the dynamics of the scalar field driving the acceler- 
ation of the early Universe. This formafism is based on a 
Langevin equation which describes the evolution of a spa- 
tiaUy averaged field (typicaUy over a Hubble patch), the 
so-cafled "coarse-grained" field. Solving this equation, 
especially when the backreaction of this coarse-grained 
field on the background geometry is taken into account, 
is notoriously known as a difficult task and various meth- 
ods have beenproposed in the literature, see for instance 
Refs. USSIi. 

Recently, a method based on a perturbative expan- 
sion in the noise was presented At second or- 
der in the noise, this method is powerful enough to en- 
sure the calculation of the probability density function 
of the coarse-grained field for arbitrary potentials. It 
was demonstrated that, in order to obtain explicit ana- 
lytical expressions, the calculation of only one quadra- 
ture is necessary. If, in addition, the volume effects are 
also determined, then only one more quadrature is re- 
quired. It turns out that these quadratures are feasible 
for a large class of inflationary models, for instance in the 
cases of the chaotic new 0, hybrid 14] and run- 
ning mass |il5| scenarios. The stochastic effects in these 
models were studied in Ref. To], where the evolution of 
the corresponding distributions was discussed in detail. 

An important question concerns the domain of validity 
of the perturbative expansion used in order to obtain the 
above results. The aim of this article is to develop a 
method to treat this question and to estimate when the 
perturbative expansion gives reliable results. It is worth 
noticing that so far (and this is also valid for the other 



approaches used in the literature to solve the Langevin 
equation) this issue has never been addressed elsewhere. 
In general, the approximate expression for the probability 
function is derived without worrying about its accuracy. 
It will be shown that the method of Ref. 0| gives, most 
of the time, a better approximation than it could naively 
be guessed on general grounds and only breaks down at 
the very end of inflation. 

Our article is organized as follows. In the next section. 
Sec. we briefly recall the main results and equations 
obtained in Ref. [T^- Then, in Sec. IIIII we present our 
method to study the accuracy of the perturbative expan- 
sion and apply it to the inflationary models discussed in 
Ref. T^, to wit, chaotic, new, hybrid and running mass 
scenarios. Finally, in Sec. lIVL we present our conclusions. 



II. SOLVING THE LANGEVIN EQUATION 

In stochastic inflation, one is interested in the behav- 
ior of a coarse-grained fleld Lp obtained after taking the 
spatial average of the original inflaton fleld over a volume 
the size of which is of the order of a Hubble patch. The 
coarse-grained field obeys a Langevin equation that can 
be written as 
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where V is the inflaton potential and ^ a white noise 
defined such that its correlation function simply reads 
(C(t)C(i')> = ^{t - t'), Hz) being the Dirac distribu- 
tion. The backreaction can be seen in the fact that 
the Hubble parameter H in Eq. ^ depends on the 
coarse-grained field if via the slow-roll Friedmann equa- 
tion ^ kV{(p)/3, where k = 87r/m^^. 

The method proposed in Ref. (see also Ref. 
for earlier attempts) consists in expanding the coarse- 
grained field in powers of the noise according to 
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where ipd is the classical solution, i.e. the one obtained 
when the noise is "switched off" in the Langevin equa- 
tion. The quantities Skpi and 5lp2 are respectively first 
and second order in the noise. They obey the equations 
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where a prime denotes a derivative with respect to the 
field. These equations can easily be solved since they are 
(by definition) linear. The expansion could of course be 
pushed to higher orders if necessary. 

In Ref. it was demonstrated that this formalism 
allows us to calculate the probability density function 
of the coarse-grained field in a single Hubble patch. At 
second order in the noise, it is given by 
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where the variance (SipY) and the mean {5tf2) appearing 
in the above expression can respectively be written as 
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This formalism also permits the calculation of volume 
effects. If, instead of considering the distribution of the 
field in a single domain, wc want to have access to its 
spatial distribution, one must weigh the single domain 
distribution by the volume of each Hubble patch. This 
leads to the definition 

(<5((p-(p[e])e3/ drH(^K])\ 



> 



Then, it was shown in Ref. that, at second order in 
the noise, Pv takes the form 
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where {Stpi') and (ip) = ipci + (^¥'2) are still given by 
Eqs. © and JT)). The term I^J describing the correction 
to the mean value due to volume effects can be written 
as 
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Therefore, as already mentioned, estimating the volume 
effects merely requires the calculation of one additional 
quadrature. 

In Ref. T(?. the results briefly described above have 
been applied to various concrete inflationary models. In 
particular, the potential 



V{ip) = 



(11) 



where a — 0,1 and b — ±1, has been considered. The 
case a = 0, 6 = 1 corresponds to large field (LF) models 
(or "chaotic inflation") [ij, a = 1,6 = —1 to small field 
(SF) models (as "new inflation") [iJl and a = 1, 6 = 1 
to hybrid inflation The scale M is fixed by the 

Wilkinson Microwave Anisotropy Probe (WMAP) nor- 
malization. The case of running mass (RM) inflation [13, 
namely 



V{ip) = 



(12) 



was also treated. In the expression of the potential, 
Afpi = TOpj/VSTT and the quantities c (which can be 
positive or negative) and po are free parameters. Run- 
ning mass inflation can be realized in four classical ver- 
sions and stochastic effects have been studied in Ref. ^3 
for the first (c > 0, ipc\ < Po) and the second (c > 0, 
fci > fo) scenarios (RMl and RM2) . 

For the models described above, the behavior of Pc and 
Pv have been investigated in details in Ref. [l3|, see in 
particular Figs. 2 and 3. As mentioned in the introduc- 
tion, the issue that we now address is the reliability of 
the method of approximation used in order to establish 
these results. 



III. RELIABILITY OF THE EXPANSION 

An important question is the determination of the in- 
terval in which the perturbed solution of the Langevin 
equation that we have obtained, (/Jci + i^iySi + 6p2, remains 
a good approximation of the exact one. Indeed, initially, 
the perturbed solution is "by definition" a good approxi- 
mation since we have i5v'i('^in) = (5(^2 (Vin) = 0. Then, as 
the field evolves from pin, we expect Spi and Sp2 to grow 
and the approximation to break down at some value of 
</'ci 7^ '/'in- j4 priori, the criterion of validity is simply 
6p2 < Spi < Pel- But things can be more complicated. 
For instance, let us assume that the classical field is ini- 
tially very small, as is the case for new inflation. Then, 
Spi/pci becomes large very quickly (because pci is very 
small), apparently signaling a breakdown of the approx- 
imation. However, it is clear that this could just be an 
artifact of the criterion used which, somehow, would be 
too naive. To illustrate this last point, let us consider the 
following simple example. Suppose that we want to cal- 
culate /((/9ci + A(/3), where / is a given function that we do 
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not need to specify explicitly. Taylor expanding this ex- 
pression leads to f{ipc\ + Ai/j) ~ fifd) + f'{'fci)^'P and, 
in general, this expression gives a good approximation 
provided that Aip <^ ipd- However, if the derivatives of 
/ are very small around fci, then the approximation can 
be good even if Aip is much larger than could naively be 
expected. We will see that, in the case of the perturba- 
tive expansion of the Langevin equation, we are exactly 
in this situation. The deep reason for that is the fact 
that the role of the derivatives of / is now played by the 
derivatives of the Hubble parameter. Since these ones are 
necessarily small as long as the slow-roll approximation 
is satisfied, one can expect the previous phenomenon to 
happen. Therefore, it is important not to underestimate 
the reliability of the perturbative expansion and to study 
this issue carefully. 

In the following, we address this question from a 
slightly different point of view, focusing on the Langevin 
equation itself rather than on its exact solution which is, 
of course, unknown. Our goal is to find a criterion that 
controls when the perturbed equation that we are able to 
solve is a good approximation of the exact one. This is a 
simpler task since we now compare known "objects" . At 
this point, one can even dare an analogy. The situation 
under consideration is indeed similar to what is done with 
the slow-roll approximation, for instance for the Klein- 
Gordon equation. In this case, one does not compare the 
exact solution (which is, most of the time, unknown as 
well) to the slow-roll one. One rather studies how small 
the term that we neglect in the exact equation (ip) is in 
comparison with the term that we keep (Hip), i.e. we 
study the magnitude of the slow-roll parameter (/3/ {Hip). 
The spirit of the method that we use below is along the 
same line. Finally, before embarking on the discussion 
of the reliability of the approximation used here, let us 
stress again that, so far and despite its importance, this 
question has not been given a satisfactory answer in the 
literature on the subject. 

In order to determine the accuracy of the expansion, 
we will make use of the Lagrange remainder theorem pTj 
for the error in a Taylor expansion. This theorem states 
that any function f{ip) around some value ipd can be 
written as 



equation up to second order in the noise. This gives 



f{ipci + Aip) 



fc=0 



(A^)", (13) 



for some value of the parameter 6 between and 1. Let 
us emphasize that this expression is exact and does not 
assume anything on Aip, in particular does not assume 
Aip < 1. 

The next step is to apply this theorem to the Langevin 
equation, ip + 2H' /n — i7^/^^/(27r), more precisely to the 
function H' and H^^^ in the left and right hand sides 
respectively. In our case, we take n — 3 since we have 
considered the perturbative expansion of the Langevin 
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where we have used the classical equation of motion and 
where, according to Eq. H13|l . we have 



(15) 
(16) 



We stress again that, despite its resemblance with 
Eqs. (0 and l@J, Eq. H14|l is an exact equation deter- 
mining Aip (hence the exact stochastic field ipd + Aip), 
as long as some values of the two parameters 0l and 0r 
are suitably chosen between and 1. At this stage, this 
is just a complicated way to re-write the exact Langevin 
equation 

The main idea is now to assume that the truncated ex- 
pansion is reliable for values of Aip = 5ipi+Sip2 such that 
L2 and i?2 are small in comparison with the other terms 
appearing in Eq. H14|l . Indeed, if this is the case, then the 
approximated Eqs. and Q become indistinguishable 
from the exact one H14|l . More precisely, for each value of 
ipd, we have to find the limiting values Aif^i^^ipd) < 
and A(y3,„ax('/'ci) > such that L2 and R2 are small in 
comparison with the other terms in the equation of mo- 
tion. Then, the validity of the perturbative treatment 
will be guaranteed as long as 

-|A(pinin(</5cl)| < Sfl + < A(/?inax(<^cl) , (17) 

or, in other words, as long as we have ip € 
[ipci-\Aip^in{ipd)\,'Pd + Aip^i,^{ipci)]. In practice, 
since we are dealing with stochastic quantities, instead 
of Sipi + 5ip2, we will apply our criterion to the quantity 
\/ (Sipi) + {Sip2), {Sip2) being evaluated with or without 
the volume effects. 

However, to explicitly derive A(^niax(95ci) and 
Aiy9niin(¥'ci), Eq- H14|l cannot be used directly because we 
do not know the values of 6'l and 6^. In fact, it is suffi- 
cient to take the maximum of the absolute value of the 
Lagrange remainders (for ^l.r G [0, 1]) in order to get an 
upper bound on the error. Therefore, the approximation 
is reliable, i.e. L2 and R2 are negligible, when the two 
following conditions 



max 
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hold, while it breaks down when (for fixed values of ipc\) 
at least one of the two is violated. The limiting values 
A</3inax(<(2ci) and A(^i„in((^ci) are then determined by re- 
quiring that the two above inequalities become equalities. 
Since we have two equations and each of them involves 
absolute values, this gives two positive and two negative 
solutions, the actual value of A(^,„ax(¥'ci) and A(^,„in((pci) 
clearly being the one leading to the tightest constraint. 

Having determined Acp^s^xi'fici) and A(pmin{(pci) with 
the above procedure, one must also take into account the 
fact that we are dealing with stochastic quantities. In 
this respect, the validity of the perturbative treatment 
will be guaranteed as long as the probability of finding 
-|A((!?inm(i,^'ci)| < Sipi +S(p2 < Av3i„ax(</'ci), IS Sufficiently 
close to 1. In terms of probability, this means that one 
requires 
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where we have considered Pc as the probability density 
function. In the case where the volume effects are taken 
into account, Pv should be used instead. 

Finally, there is yet another constraint coming from 
the fact that, in general, 5ipi +5ip2, is a good approxima- 
tion only if 5(p2 ^ <5v?i • This is necessary if we want to 
"separate" Eq. H14|) into two equations, one for 5ifi and 
one for 5ip2- 

Let us now see how the previous considerations work 
in practice for the chaotic infiation potential V{ip) = 
m^(/3^/2. In this particular case, dipi + 5ip2 is an ex- 
act solution of the approximated second order equa- 
tion since H'" = and the constraint S(p2 ^ Sipi 
does not apply. In addition, we also have L2 — 
and, as a consequence, the limiting values Aip^i^ and 
Aiy9max are found only from the constraint H19|) involv- 
ing R2. Using the slow-roll equations of motion, one has 
(i?3/2)'(a.) = 3(«;/6)^/V3/2a:i/2/2 and (H^/^)" {x) = 

The next step is to evaluate 



3(K/6f/V3/2x-i/V4. 



the maximum of this last function in the interval x e 
[ipci, ifci + Aip]. Let us start with the upper bound. Since 



> one has max 
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Then one can solve for Aip„ 
arrives at 
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I.e. 



0. 



Applying Eq. ((T^ . one 
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from which one obtains A(^niax = 4(/Jc1- 

Let us now consider the lower bound A(^min < 0. 
The maximum of the function is now given 



by max 



-1/2 



CX (<^cl - \Alfni 

Therefore, in this case, solving the corresponding equal- 
ity (|19|) requires to solve a second order algebraic equa- 



tion in |Aiy9rnin|, namely 

1 3 /K\3/4 



4 V6 



m3/2(<^,i- |A(p„,i„|) |A</j,„i„|2 
^(^)'''-^/VcY>^n.i„|, (22) 



and the result reads |Ac/j„iin| — (—8 ± 4\/5)<^ci- Gath- 
ering the two limits obtained before, one finds that the 
reliability interval is given by 



-4(\/5 - 2)(/7ci < 5ifi + 5(p2 < 4(pc 



(23) 



This interval is represented in Fig. (left panel) by the 
uniformly colored (grey) region. One clearly sees that 
this region is limited by two straight lines as calculated 
above. As inflation proceeds, the allowed re gion s hrinks. 
This has to be compared with {S1P2) ± \/ repre- 
sented by the hatched (blue) area. The lower border 
line of the hatched (blue) region is {8(^2) — \/ (S^i) while 
the upper border line is {6lp2) + \/ i^'Pi)- The two lines 
meet at the beginning of inflation where they vanish since 
(^S(pl{(pin)) = {Sip2{<Pin)} = 0. As long as the hatched re- 
gion lies within the uniformly colored one, the approx- 
imation is reliable. When this is no longer the case, 
the approximation breaks down. In Fig. ^ this is sig- 
naled by the vertical dotted-dashed line and occurs for 
(fci ^ 2 X lO^mp,. It is clear that the second order ap- 
proximation is good until one approaches the end of slow- 
roll inflation. The right panel of the same figure shows 
the probability of finding dipi + Sip2 between Aip^i^ and 
Aiy9max, computed according to Eqs. (^0)) and (j^Sl, and 
confirms the previous conclusion. 

When volume effects are considered, the situation be- 
comes conceptually more complicated but the same ideas 
can be utilized to check the accuracy of the volume- 
weighted distribution. In particular, one shoul d now 
compare the region limited by {6lp2) + 'il^J ± \/ {StpD 
with the reliability region. In the case of large field mod- 
els, however, we do not plot this region because the vol- 
ume effects are so important that the corresponding re- 
gion would be outside the figure. This will be done for 
the other models, see below. 

One can also use the criterion of Eq. H2()(l but this time, 
as already mentioned above, the probability should be 
evaluated with the distribution Pvi^) rather than with 
Pc{^). The corresponding results can strongly differ 
since the field realizations having higher potential en- 
ergy (i.e. with a faster expansion rate) will be favored. 
In particular, if their expansion rate is sufficiently large, 
this can give a high statistical significance to realizations 
outside the reliability range having a very low signifi- 
cance according to the original distribution. In this situ- 
ation, when the difference between the two distributions 
is very important, the form of Pv{^) obtained from the 
perturbed solution cannot be trusted although Pc(</') is 
reliable. This means that, in most domains, the statisti- 
cal properties of the field are correctly described by Pc{f) 
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FIG. 1: Accuracy of the second order approximation for the large field (chaotic) model V{(p) oc rri^ifp' with an initial condition 
corresponding to V(v5in) = mpj/2. The mass m is chosen so that the WMAP normalization is reproduced (its value is related 
to the value of M used before). All quantities are plotted versus ^ci and, therefore, inflation proceeds from the right (large 
field values) to the left (small field values). On the left panel, the allowed interval is represented by the uniformly colored 
(grey) region which is delimited by A(/Pmin and Ai^max obtained with Eq. (1231 . The hatched (blue) region represents the region 
delimited by the two lines (<5(p2) ± \/ {Sifil), (5</32) being evaluated without the volume effects. The vertical dotted-dashed (blue) 
line signals the value of <pci at which the approximation breaks down. On the right panel, the probability of finding ip in the 
reliability range computed with Pc (solid blue line) and with (red dashed line) is displayed. Clearly, the single-domain 
probability probability starts decreasing approximatively at the value of i/Jci where Sifi + 5ip2 is no longer in the reliability 
interval. On the other hand, the volume effects corrections are very large and the volume weighted distribution is not trustable. 
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FIG. 2: Accuracy of the second order approximation for the small field (new inflation) potential V{if>) oc 1 — {(p/fi)'^ with the 
initial condition (fin — 10~^fi. In this case, inflation proceeds from the left (small fleld values) to the right ("large" field values). 
On the left panel, the uniformly colored (grey) region represents the interval where the approximation is trustable. The limiting 
values are now obtained from the condition on the remainder R2 but also from the one coming from the remainder L2. The 
actual A(^min(v5ci) and Ai/5max(v5ci), which delimit the confidence region, must be the smallest ones (in absolute value). The 
hatched (blue) region with positive slope lines is delimited by {5ip2) ± yJSipf) where {Sip2) is computed without the volume 
effect. On the other hand, The hatched (red) region with negative slope lines is also delimited by (^v'a) ± \/ (Sfi) but, this 
time, with {Stfi2) computed with the volume effects. The vertical dotted-dashed (blue) line indicates when the approximation 
without the volume effects breaks down while vertical dotted (red) line signals when the approximation with the volume effects 
becomes untrustworthy. On the right panel, the probability of finding Sifi + Stfi2 in the reliability range is displayed. It is clear 
that both Pc (solid blue line) and Pv (dashed red line) yield a probability close to one for a large part of the inflationary phase. 



but that the Universe is mainly made up of very big do- 
mains where Pd^fi) cannot be trusted. This is exactly 
what happens for LF models where the volume-weighted 
probability of finding the field in the confidence range 
(dashed red line, right panel in Fig. ^ is basically vanish- 
ing while the single-domain one (solid blue line) is large. 
Therefore, in this case, the perturbative method does not 



allow us to reliably compute the volume-weighted distri- 
bution. 

We have also performed the same study for new infla- 
tion and the results are displayed in Fig. |21 Now, the 
reliability interval, still given by the uniformly colored 
(grey) region, is no longer determined only from R2 but 
also from L2 because, contrary to the LF case, L2 does 
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FIG. 3: Accuracy of the second order approximation for the two running mass models RMl and RM2 with initial conditions 
Vi-a/fo — 1 — 1.5 X 10~^ and (fin/^o — 1 + 10"'^ respectively. Inflation proceeds from the right to the left in the RMl model 
(upper panels) and from the left to the right in the RM2 model (lower panels). On the left panels, the limiting values signaling 
the break down of the conditions given by Eqs. H15|l and H16|l are represented and compared with {5ip2) ± \/ {S(pi)- Conventions 
are the same as the ones used in Figs. and |5| As usual, the actual values of A<^min and A<^max must be the ones with 
the smallest absolute values. Cusps in the curves are a consequence of taking the maximum absolute value of the Lagrange 
remainders and appear when one discontinuously change the value of from to 1 (or the opposite). On the right panels, the 
probability of flnding ip in the reliability range is displayed. For both models the reliability of the solution does not dramatically 
change (and for the RM2 model is even slightly enhanced) when volume weighting is considered 



not vanish for SF. Therefore, in this situation, one has 
to determine Aipj^a^ and A(/3i„in from Eqs. ((TH|l and ifT^ 
and not only from Eq. H19|l . as was the case for the LF 
models. The form of A(/3niax and Atpmin as a function 
of ipci is also more complicated and is no longer given 
by straight lines. It must be computed numerically. As 
can be seen in Fig. 13 the allowed region increases at the 
beginning of inflation, reaches a maximum extension and 
shrinks as the end o f inflation is approached. The re- 
gion {S1P2) ± \/ i^fi) without the volume effects is given 
by the hatched (blue) region, the lines having a positive 
slope. The region {S(p2) ± \/ {5^p\) with the volume ef- 
fects taken into account is represented by the hatched 
(red) region, the lines having a negative slope. When 
these regions are within the uniformly colored (grey) re- 
gion, the approximation is reliable. As before, the border 
lines {5lp2) ± \/ {S^pf} meet at the beginning of inflation 
where they vanish. The break down of the approxima- 
tion, without the volume effect, is signaled by the vertical 
dotted-dashed (blue) hue and, with the volume effects. 



by the dotted (red) line. One notices that the results for 
SF are basically similar to those that have been obtained 
in the chaotic model case, namely the approximation re- 
mains reliable until the very end of the inflationary phase. 
However, in the case of new inflation, one clearly sees the 
importance of using a carefully defined criterion to esti- 
mate the reliability of the approximation. As already 
mentioned, the naive criterion Sipi <C (pd would have in- 
dicated that the approximation becomes untrustworthy 
very quickly after the beginning of inflation since we have 
initially ipd/fJ. ^ 1. We see in Fig. |21 that, on the con- 
trary, the approximation is good during a large part of 
inflation. Finally, the probability of being in the relia- 
bility region, computed with Pc (solid blue line) or 
(dashed red fine) is also displayed in Fig. |3 (right panel). 
In the case of SF models, the volume effects are less im- 
portant and, as a consequence, the two probabilities are 
similar. The fact that the break down of the approxima- 
tion occurs at the end of inflation only, at </?ci/a* ~ 0.65 
for the particular example studied here, is confirmed. 
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A similar analysis can also be done for the running 
mass potential. The results for the two models under 
consideration are displayed in Fig. 13 All the conventions 
concerning the allowed regions, volume effects etc . . . are 
the same as before. In the case of the RMl model (upper 
panels), Eq. H18(l giving the constraint on L2 can be an- 
alytically solved to the lowest order in c((y9o/MpJ^. This 
gives 

< Av3 < 3(^ci . 24 

D 

It turns out that this constraint is the dominating one at 
late times as can be checked in Fig. O (left panel): the 
shape of the allowed region near the end of inflation is 
delimited by two straight lines. The constraint on R2 
has been solved numerically. The figure demonstrates, 
and this is also confirmed in the right panel, that the 
perturbative solution for RMl is very good during almost 
all the inflationary phase and breaks down only at the 
very end of inflation. Similar conclusions hold for the 
RM2 model, see the two bottom panels. 

Finally, one notices that the two probabilities [i.e. the 
ones obtained with /'c(v) ^^^d Pv(¥')], and contrary to the 
LF case, do not dramatically differ from each other. In 
the case of the RM2 model, the reliability of the volume- 
weighted description can even be larger than the single- 
point one. 

IV. CONCLUSIONS 

In this section, we briefly summarize the new results 
obtained in this article. The main goal of the paper was 



to present a new method aimed at estimating the preci- 
sion of the perturbative expansion studied in Ref. |lClj| . 
This method is based on the use of the Lagrange re- 
mainder. After having discussed the general features of 
this new approach, we have applied it to the inflation- 
ary models studied in Ref. ^3 ■ We have proven that the 
approximate probability density functions derived in this 
reference are, in general, a very good approximations to 
the actual ones except, as expected, at the end of infla- 
tion. This conclusion holds even if the volume effects are 
taken into account except in the case of the large field 
models. We conclude that the perturbative expansion 
of the inflationary Langevin equation together with the 
method presented here, besides being the only available 
method with a built-in measure of its domain of valid- 
ity, form a robust formalism to efficiently compute the 
stochastic effects during inflation. 
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